Out-of-equilibrium mean-field dynamics of a model for wave-particle interaction 
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The out-of-equilibrium mean-field dynamics of a model for wave-particle interaction is investi- 
gated. Such a model can be regarded as a general formulation for all those applications where 
the complex interplay between particles and fields is known to be central, e.g., electrostatic insta- 
bilities in plasma physics, particle acceleration and free-electron lasers (FELs). The latter case is 
here assumed as a paradigmatic example. A transition separating different macroscopic regimes is 
numerically identified and interpreted by making use of the so-called violent relaxation theory. In 
the context of free-electron lasers, such a theory is showed to be effective in predicting the satu- 
rated regime for energies below the transition. The transition is explained as a dynamical switch 
between two metastable regimes, and is related to the properties of a stationary point of an entropic 
functional. 



I. INTRODUCTION 

Mean-field models have been widely studied as 
paradigmatic representatives of the important class of 
systems subject to long-range coupling. In the simplest 
scenario, particles are made to interact in one dimen- 
sion, subject to a varying field which is self-consistently 
sensitive to the individual trajectories. A global network 
of connections is hence driving the dynamics of every 
constituting element, as it certainly happens for more re- 
alistic settings where e.g. gravity or unscreened Coulomb 
interactions are at play. 

The dynamics of mean-field models displays intrigu- 
ing features. Particles may be trapped in intermediate 
(out-of-equilibrium) states, whose duration diverges with 
the number of constitutive elements, and which substan- 
tially differ from the corresponding thermodynamic equi- 
librium configuration. These metastable states are often 
termed in the literature Quasi-Stationary States, here- 
after QSS, and bear an extraordinary conceptual impor- 
tance as they potentially corresponds to the solely exper- 
imentally accessible regimes, for a wide range of applica- 
tions ranging from celestial mechanics to plasma physics. 

Interestingly, the evolution of the QSS is intimately 
governed by the discreteness of the medium being investi- 
gated. More specifically, QSS are believed to correspond 
to stationary stable solution of a Vlasov model, invoked 
as the continuous analogue of the discrete A^-particles 
dynamics. 

Within this context, the Hamiltonian Mean Field 
model (HMF) [H has often been referred to as the bench- 



mark model for elaborating onto the QSS emergence. 
This is a one-dimensional Hamiltonian describing the 
evolution of rotors coupled via a mean-field cosinus- 
like potential. The QSS in the HMF setting have been 
explained by resorting to a maximum entropy principle, 
pioneered by Lynden-Bell in astrophysical context, and 
fidly justified from first principles |^. Here, the sup- 
posedly relevant Vlasov picture enters the description as 
a Fermionic contribution to an entropy functional. Be- 
sides, out-of-equilibrium phase transitions are predicted 
to occur, separating between distinct macroscopic QSSs. 

The Lynden-Bell protocol, also termed violent relax- 
ation theory, was also argued to apply to other mean-field 
models and, indeed, it showed effective in predicting the 
saturate intensity of a free-electron laser (FEL). FELs are 
lasing devices consisting of a relativistic beam of charged 
particles, interacting with a co-propagating electromag- 
netic wave. The interaction is assisted by the static and 
periodic magnetic field generated by an undulator. FELs 
admit a mean-field description in term of the so-called 
Colson-Bonifacio modelj^Q, which captures the essence 
of the collective wave-particle dynamics. However, no de- 
tailed study has been carried out for the FEL case aiming 
at unravelling the possible existence of out-of-equilibrium 
transition of the type mentioned above. Are these transi- 
tions ubiquitous in mean field dynamics and, in this case, 
can we provide a consistent intrepretative framework for 
their emergence? This paper is dedicated to answering 
such questions, and makes reference to the specific FEL 
setting. We also stress that the Colson-Bonifacio model 
of FEL dynamics can be regarded as a general formula- 
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tion for all those applications where the complex inter- 
play between particles and fields is well known to be cen- 
tral, e.g. electrostatic instabilities in plasma physicsQ- 
As a closing remark, we notice that, on the practical im- 
plication side, by disposing of reliable predictive tools on 
the system evolution, one can aim at guiding the system 
towards different experimental regimes. 



II. THE FEL MODEL 

The Colson-Bonifacio model for the FEL dynamics de- 
scribes the coupled evolution of the electrons with a co- 
propagating wave. The equations read: 



dpj_ 

dz 

dA 



A*e- 
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where 9j stands for the particle phase with respect to 
that of the optical wave, Pj being its conjugate normal- 
ized momentum. The complex quantity A = + iAy 
represents the transverse field and N the number of elec- 
trons composing the electron bunch [27j . 

Here z labels the longitudinal position along the un- 
dulator, and it effectively plays the role of time. The 
intensity of the laser field is / = -I- Ay. As it can 
be seen from the last of Eqs. (tlj, the bunching term. 



is the source of wave amplification. The 



bunching quantifies the degree of localization of the elec- 
trons in the generalized space of their associated phases. 
The above discrete system of equations admits a Hamil- 
tonian formulation to which we shall make reference as 
to the N-body model. In the N ^ oo limit, system ^ 
converges to the following Vlasov-wave set of equations 
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Eqs. ([2]) can be simulated numerically, thus allowing 
us to monitor the evolution of the phase space distri- 
bution function f{9,p) along the z axis. In our imple- 
mentation, we adopt the semi-Lagrangian method 0, 
associated with a cubic spline interpolation jsj. Re- 
sults of the numerical integration are also checked ver- 
sus A^-body simulations and shown to return a perfect 
matching on relatively short time scale, for large enough 
values of N. On longer times, finite-TV corrections do 
matter. The discrete system is in turn sensitive to in- 
trinsic granularity effects, stemming from the intimate 
finiteness of the simulated medium, and progressively mi- 
grate from the Vlasov state towards the deputed equilib- 
rium configuration. When increasing its size, the system 



spends progressively more time in the Vlasov-like, out- 
of-equilibrium regime. Formally, in the iV — >■ oo limit, it 
never reaches equilibrium, being permanently trapped in 
the QSS. 

As previously anticipated, our study is hence ulti- 
mately concerned with the emergence of QSSs, in a con- 
text where particles and waves evolve self-consistently. 
We shall be in particular interested in elucidating the oc- 
currence of out-of-cquilibrium phase transitions via dedi- 
cated numerical simulations, and substantiate our claims 
analytically. In doing so, we will virtually extend the 
conclusion of Q to a broad spectrum of potentially rel- 
evant applications, beyond the specific case under in- 
spection. Among others, it is again worth mentioning 
plasma physics: A formulation equivalent to model ^ is 
in fact often invoked, when studying the collective effects 
of beam-plasma dynamics jsj . 



III. ON THE INITIAL CONDITIONS AND 
THEIR SUBSEQUENT DYNAMICAL 
EVOLUTION 



Let us turn to discussing our results, as obtained via 
numerical integration of ([2]). In order to make contact 
with the investigations reported in we shall employ 
in the following a two-dimensional water-bag initial con- 
dition in phase space, which can be seen as a rough ap- 
proximation of a smooth Gaussian profile. A (rectangu- 
lar) water-bag is formally parametrized by two quanti- 
ties, namely the semi- width of the spanned interval in 
phase, A9, and its homologue in the momentum direc- 
tion, Ap. The corresponding expression for / can be cast 
in the form (see also Fig. [S] top- left): 



fiO^p) 



fo if \p\ < Ap, 
\9\<A9, 
otherwise. 



(3) 



The initial conditions can be also characterized by defin- 



ing 



e = 



I Ad 



Ap 
6 



(4) 



where &o is the initial bunching, and e the initial average 
kinetic energy per particle. Notice that we thus access 
all possible values of the bunching bo G [0, 1] by properly 
tuning A9, and all positive energies e by varying Ap. We 
here limit our discussion to the case of vanishing initial 
optical field, Iq ~ 0, the relevant parameters' space be- 
ing therefore solely bound to the plan (6o, e)- The initial 
condition here selected is hence solely controlled by these 
two parameters. In other words, we are specializing on 
a given bidimensional subset, and deliberately ignore the 
third, in principle available, direction of the reference pa- 
rameter space. Quantifying the role of such an additional 
degree of freedom ultimately amounts to investigate the 
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so-called seeded configuration [T^, [l^ and will be the 
subject of future work. 

Let us start by discussing the simplest scenario, where 
the initial beam of particles is uniformly distributed over 
[— 7r;7r]. From a physical point of view, this amounts 
to specialize to the case of Self- Amplified Spontaneous 
Emission (SASE, where bo = and no seed is 

applied externally. Such a choice was also considered by 
Barre et al. Q and Curbis et al. [l^, where the de- 
pendence of the system evolution on the energy was nu- 
merically monitored, within the iV-body discrete view- 
point. Interestingly, 69 = is a stationary solution of 
the Vlasov system: a local perturbative calculation can 
hence be straightforwardly implemented so to investigate 
its inherent stability [ll| ; the calculations are detailed in 
appendix 1X1 For e < 0.315, an instability occurs: both 
the wave intensity and the bunching factor rapidly grow, 
before relaxing towards an oscillating plateau. The aver- 
age value of / reached in the oscillating regime is called 
the saturated intensity /. This behaviour is displayed in 
Fig. [TJ where the simulations with e > 0.315 do not show 
an amplification of /. This is a well-known property, in- 
deed correctly reproduced by our numerical simulations, 
and which first signals the existence of phase transitions, 
of the type depicted in Q. To further corroborate our 
guess on the 60 = behaviour, we turn to measuring the 
saturated intensity / as function of the energy e, where 
/ stands for the mean of / after saturation. It is here 
computed during four oscillations of /. 

As shown in Fig. [21 / rapidly shrinks, when increas- 
ing the energy e, until a critical value is reached where 
a sudden transition to / ~ is observed, bearing the 
characteristic of a first order phase transition. This is a 
further point of contact with the analysis carried on in 
[il for the HMF toy model. 

Motivated by these findings, and to push the analogy 
with the HMF setting, we consider bunched initial dis- 
tributions. From a physical point of view, this choice 
is relevant to the case of FELs working in the so-called 
harmonic generation regime [l3i[3- Particles' positions 
are here initially assigned so to uniformly span a limited 
portion of the allowed support, symmetric with respect 
to the origin, controlling the associated bunching via Eq. 
dH). The inhomogeneous (60 7^ 0) distribution in phase 
space is by nature non-stationary. Vlasov dynamics can 
however smooth it to a homogeneous, (6 = 0, / = 0) pos- 
sibly non water-bag, distribution or evolve to a bunched 
situation. The saturated mean-field average intensity / 
vs. the energy parameter e is represented in Fig. [21 show- 
ing the newly collected data for different values of 60 > 
to the reference profile relative to 60 = 0. In all cases the 
intensity is shown to decrease, as the energy increases. 
Importantly, for small values of 60, an abrupt transition is 
observed, which can be naively interpreted as of the first- 
order type. For larger values of ho, the observed transi- 
tion becomes smoother, such as for a second-order one. 
A substantially identical scenario holds for the bunch- 
ing, which evolve towards an asymptotic plateau 6, also 













£ = 0,05 






1.2 




- - c = n,i 






1,0 




E = 0,2 

- - t = 0-,S2 
■ ■ E = 0,1 




f 


0.8 










0.6 










0.4 










0.2 








0,0 









20 40 60 80 100 

Z 



FIG. 1: (Color online) As a result of the numerical integration 
of Eqs. ([2]), the evolution of / as a function of z is reported for 
different choices of the energy e, (e = 0.05,0.1,0.2,0.3,0.4); 
bo = 0. Symbols pinpoint the position of the first peak in 
the intensity time series, thus returning an indication on the 
saturation time. Notice that for e > Ec = 0.315, the peak 
is found for / <C 1 : the corresponding initial conditions are 
hence stable, and no instability develops. 




FIG. 2: (Color online) Saturated intensity, /, 
vs. e different choices of the initial bunching 

60 = 0.0, 0.05, 0.20, 0.50 and 0.90 

sensitive to the e and 5o parameters, see Fig. [3l This 
scenario points towards a unifying picture on the emer- 
gence of out-of-equilibrium phase transitions within the 
considered class of mean-field Hamiltonian model. As 
previously anticipated, the Lynden-Bell theory of violent 
relaxation was successfully applied to the HMF problem, 
allowing one to gain a comprehensive understanding on 
the out-of-equilibrium phase transition issue, including a 
rather accurate characterization of the associated transi- 
tion order. In the following section we set down to apply 
the Lynden-Bell argument to the present case, bench- 
marking the theory to numerical experiments. 
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FIG. 3: (Color online) Same as Fig. [5J for the saturated 
bunching, b. 



Before ending this section, we briefly discuss the 
phase-space structures resulting from the Vlasov-based 
simulations. Two phase-space portraits are enclosed 
in Fig. 01 and refer to different values of the energy 
e, respectively below (upper panel) and above (lower 
panel) the critical transition energy relative to the se- 
lected (fixed) 60 amount. When the system evolves to- 
wards a state at / 7^ 0, then f{9,p) shows a large res- 
onance. At variance, in the opposite regime, a hole- 
rcsonancc dipole structure is observed, see also [Ifll- This 
observation seems to suggest that the out-of-equilibrium 
phase transition materializes via a bifurcation of invari- 
ant structures, an observation that was recently made for 
the HMF model • Results of further investigations on 
this specific topic will be presented in a separate contri- 
bution [11. 
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FIG. 4: (Color online) Phase space density f{6,p) for bo 
0.05 and e = 0.2 (top) and 0.4 (bottom). 



IV. ON THE VIOLENT RELAXATION THEORY 

In his work on self-gravitating systems, Lynden-Bell 
suggested that the collisionless dynamics governed by 
the Vlasov equation tends to maximize a Fermionic en- 
tropy. The latter is obtained from the classical defini- 
tion, where the counting of the microscopic states, com- 
patible with a given macroscopic configuration, results 
from a combinatorial calculation, and is sensitive to the 
underlying Vlasov dynamics. The method was success- 
fully employed in the study of the HMF model [i,[ll,[23| 
and also applied to predict the quasi-stationary ampli- 
tude of the FEL wave 0, [l^. In these works, however, 
the analysis just focused on the unstable regime (/ 7^ 0): 
no attempt was in fact made to reconcile it, with the 
high energy homogeneous state, via the phenomenon of 
out-of-equilibrium phase transitions. More recently, Ya- 
maguchi [2l| used the Lynden-Bell approach to predict 
the core of the gravitational sheet model, demonstrating 
its adequacy within the field for which it was originally 
conceived. In the following we shall review the main steps 
of the derivation of the violent relaxation theory, applied 
to the FEL setting. Starting from a water-bag, the en- 
tropy to be maximized can be cast in the form [l^l 

= - fdpde [1 In ^ + f 1 - ^) In f 1 - ^ 
J L/o Jo V JO / \ JO, 

.(5) 

where /o is specified in (jS]) and / is the coar se g rained 
distribution function. Following Barrc et al. (6|,|23|. max- 
imizing the functional ([5]), results in the following set of 
equations 
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and F2{y) = 
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Here /3 and x are (rcscalcd) 



Lagrange multipliers and ultimately stem from the con- 
servation of mass, momentum and energy. A, /3 and 
X are calculated by solving Eqs. ([6|) numerically via a 
Newton-Raphson method. The resulting (real) value of 
A is expected to return an estimate of the laser intensity 
at (Vlasov) saturation, I = A^, while f{e,p) is: 



f{0,p) = fo 



1 



(7) 



For A = (namely / 
(j6]) reduces to: 



1 + X e/5(P^/2+2^sin6l-|-A2p+Ai/2) ' 

0) the optimization problem 
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An X value exists which solves the above equation for 
any choice of A0. The homogeneous state is a stationary 
solution of the Lynden-Bell entropy, and so a potentially 
attractive state of the Vlasov dynamics. Additional in- 
homogeneous solutions {A ^ 0, or, equivalently, / 7^ 0) 
might however emerge from investigating the full system 
(pl . The homogeneous and inhomogeneous solutions will 
be referred to as to LBO and LBA, respectively (sec Fig. 
[5]) . The forthcoming discussion will focus on how to dis- 
criminate between the two, and eventually predict the 
asymptotic fate of the svstem. 
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FIG. 5: (Color online) Examples of f{d,p): (top to bottom, 
left to right) the waterbag, the LBO solution, the LBA solution 
and the velocity distribution function integrated over space 
for the three solutions. 



V. INTERPRETING THE 
OUT-OF-EQUILIBRIUM TRANSITION AS A 
DYNAMICAL SWITCH BETWEEN LBA AND 
LBO 

Let us focus first on the LBA solution. In Fig. [51 
we report the value of / as it follows from Eqs. ([6|), for 
different choices of the energy and initial bunching. The 
predicted intensity is shown to decrease, when the energy 
gets larger, but no transition is observed, in contradic- 
tion with the results of our numerical simulations. As 
previously stressed, the homogeneous LBO state is also 
solution of the optimization problem 1^ and could in 
principle prevail over the former. To shed light on this 
issue, we calculated the entropy values Sa and 6*0, as- 
sociated to LBA and LBO, respectively. Results of the 
computations are shown in Fig. [3 where the depen- 
dence on the energy e is monitored for various choices of 
bo- Surprisingly, and at odd with what happens for the 
HMF model Q, Sa is always larger than Sq. The two 
curves do not cross each other and the LBA configura- 
tion is entropically favoured. Let us note that for higher 



values of e, / decreases and the two solutions get close 
to each other, also from the point of view of the entropy, 
without crossing. 
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FIG. 6: (Color online) The saturated intensity / for the (in- 
homogeneous) LBA solution of system ((Sjl plotted as function 
of the energy e. Different curves refer to distinct value of the 
initial bunching 60 (see legend). 

The observed transition could possibly stem from a 
purely dynamical mechanism, and this would justify the 
discrepancy between the simulation output and the sta- 
tistical prediction. More specifically, we here argue that, 
depending on the selected initial conditions, the system 
explores a local basin of attraction and struggles to find 
its way to the deputed, global maximum of the entropy. 
To clarify this point, we focus on a single numerical sim- 
ulation, assuming the system to be initialized in a LBO 
state. The dynamics can progressively take the system 
towards the LBA configuration, respecting the maximiza- 
tion of the Lynden-Bell entropy. The opposite is not pos- 
sible, and a simulation started in the LBA state will cer- 
tainly not evolve to the LBO. However, dynamical effects 
might be also at play and interfere with ideal situation 
here schematized, by virtually blocking the system in the 
neighborhood of an initially assigned LBO conformation. 
Is this the correct scenario? And how can one explain the 
observed transition when starting from waterbag initial 
conditions? These issues are addressed in the following, 
where the stability of LBO and LBA is investigated via 
direct Vlasov simulations. 

In Fig. [5] the dynamical evolution of the intensity / is 
depicted, for three different classes of initial conditions, 
relative to the same choice of e and 69. The LBA is indeed 
stable, no deviation from the initial configuration being 
observed as an effect of the Vlasov dynamics. Conversely, 
the LBO condition proves unstable, and the intensity con- 
verges towards an oscillating plateau. Interestingly, the 
LBO and water-bag (WB) evolutions are qualitatively 
similar, and, moreover, display the same average asymp- 
totic value for the intensity /. Even more important, the 
asymptotic value corresponds to the LBA (maximum en- 
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FIG. 7: The Lynden-Bell entropy calculated respectively for 
the LBO (plain line) and LBA (dash-dotted line) solution. 
Here 60 = 0.00, 0.50, and 0.90 (resp. panels 1,2 and 3). 



tropy) solution. 
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FIG. 8: (Color online) The intensity / as a function of time, 
for three different choices of the initial condition: WB (plain 
line), LBA (dashed line) and LBO (dotted line). All conditions 
refer to foo = 0.05,e = 0.10 

To further elucidate the analogies between LBO and 
WB, and also clarify the stability issue, we performed 
a detailed campaign of simulations, aimed at generaliz- 
ing the results of Fig. [51 Results of the investigations 
are reported in Figs. [S] and I10[ where / is represented 
versus e, for two selected 60 amounts. The LBA initial 
condition is always stable under Vlasov dynamics. It is 
in fact a global maximum of the Lynden-Bell entropy, 
which, in this respect, proves adequate to describe the 
system at hand. The observed LBO evolution is by far 
more complex. For large values of the energy, LBO is 
stable. The stability is eventually lost when reducing 
the energy parameter: A transition materializes and the 



LBO evolves towards the LBA state. In the vicinity of 
the transition, the LBO initial condition approaches an 
asymptotic configuration, which slightly differs from the 
LBA one, and possibly results from a balance between 
opposing dynamical strengths. The WB evolution mim- 
ics that of LBO; the two curves returning a pretty close 
correspondence. In practice, during a short transient, 
an initially bunched WB expands (almost) ballistically, 
the particles being essentially transported by their own 
initial velocities, so visiting the whole interval [— 7r;7r]. 
The obtained distribution can be approximated by a ho- 
mogeneous LBO state (the field has not yet developed, 
its intensity being effectively negligible), which in turn 
explains the observed correspondence. We shall however 
emphasize that this mechanism applies to relatively small 
&o amounts. This fact is testified in Fig. (TU] The agree- 
ment between LBO and WB evolution is shown to worsen, 
when compared to that of Fig. [3] This is understood as 
follows; Starting from a high degree of bunching, the in- 
duced field opposes the natural ballistic contribution, by 
further enhancing the tendency to form a coherent clump 
of particles. 

In summary, our calculation returns two stationary 
points of the Lynden-Bell entropy. The first, which we 
termed LBA, corresponds to a inhomogeneous (laser on) 
configuration and it is a global maximum of the entropy. 
The second, labelled LBO, is homogeneous (laser off). 
For sufficiently large energies, the system can be locally 
trapped in the vicinity of the LBO. This happens if the 
system is initiated close enough to a LBO state, as e.g. in 
the case of WB with moderate 60 values. For smaller e, 
the LBO loses stability and the system departs towards 
the entropically favoured LBA state. 

Having detected no additional stationary points of the 
fermionic entropy, other than LBO and LBA, we inter- 
pret LBO and LBA as a saddle point and a global maxi- 
mum, respectively. While LBO could also be in principle 
a global minimum, this hypothesis is invalidated by the 
fact that there exists (at least) a function /. compatible 
with the system dynamics, which yields a lower entropy 
value. Consider in fact ; 

\0 otherwise. ^' 

This is not a stationary point of Lynden-Bell's entropy 
but represents one of the admissible equilibria of the dy- 
namical system ([2]), in its coarse grained perspective as 
implied by the theory. The Lynden-Bell entropy associ- 
ated to Eq. ([9|) can be straightforwardly computed via 
Eq. ([5]) and it is found to be always smaller than 6*0, 
the value associated to the LBO configuration, for all 60 
and e. Based on the above, and as previously antici- 
pated, we can exclude the possibility for LBO to be a 
global minimum. Following our deductive reasoning, we 
hence suggest that LBO is instead a saddle-point and the 
observed transition is consequently interpreted to stem 
from a local modification of LBO stability properties or 
morphological characteristics (e.g. width/flatness of the 
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stability basin) . A detailed analytical characterization of 
LBO stability is at present missing and could eventually 
help clarifying the underlying scenario. 
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FIG. 9: (Color online) The saturated intensity / is plotted 
as function of the energy e, for bo = 0. WB, LBO and LBA 
initial conditions are considered, see legend. 




FIG. 10: (Color online) The saturated intensity / is plotted 
as function of the energy e, for 6o — 0.50. WB, LBO and LBA 
initial conditions are considered, see legend. 



VI. CONCLUSIONS 

Using the case of a FEL as a paradigmatic example, in 
this paper we described the out-of-equilibrium dynam- 
ics of a mean-field model for wave-particle interaction. 
Our numerical investigation moves from the Vlasov ver- 
sion of the model, which rigorously applies to the con- 
tinuous limit and is believed to constitute the correct 
interpretative framework to elaborate on the QSS pecu- 
liarities. Working within this context, and assuming a 



specific class of initial conditions, we identified a switch 
between different macroscopic regimes. Such a transition 
is ultimately controlled by the nominal energy value and 
by the initial particles bunching, in qualitative agreement 
with what was previously observed for the HMF model. 

The Lyiiden-Bell violent relaxation theory is developed 
with reference to the FEL setting to quantitatively sub- 
stantiate our findings. We numerically characterized the 
stability of the stationary points of the Fermionic en- 
tropy functional. A sudden change in their characteristic 
is related to the occurrence of the observed transition, 
supporting the adequacy of the violent-relaxation theory 
to describe the states reached by the dynamical system. 

As a final comment, and beside stressing the unifying 
picture that is here brought forward, we emphasize that 
the transition here predicted can be in principle observed 
in real devices. We regard this as a rather important 
series of experiments which could eventually result in a 
direct proof on the existence of QSS for wave-particle 
systems. 



Appendix A: Stability of the homogeneous waterbag 

This Appendix is devoted to reviewing the stability 
condition of system ^ for an initial homogeneous distri- 
bution of the waterbag type. The derivation follows from 
a straightforward linear analysis which can be found for 
instance in [5|. We shall hereafter make reference to the 
calculation detailed in 



24| . Let us start by assuming 



a general equilibrium setting where the spatial distribu- 
tion is homogeneous {A^ = Ay = 0) and / = fo{p), i-e. a 
generic function of the variable p. Then one can linearize 
around the equilibrium and eventually derive an explicit 
solution which holds for a relatively short time. To this 
end we write: 



f{0,P,t) ^.Mp) + fi{e,p,t), 

A,{t) -Ai(t) and Ay{t)^Yi(t) . (Al) 

where the quantities labeled with the index 1 stand for 
the linear perturbation. Introducing in system (2) and 
retaining the lowest order yields: 



{dg+pde)fi-2r]{XiCos0 -^Yisine) =0 (A2) 

r /'+00 

de dp /i cos 6* =0 (A3) 

K J — OO 

/TT p-\-QQ 
de dp h sine +^ =0 (A4) 
'TT J —QO 

where we introduced ri{p) = dpfo{p). The above linear 
system admits a solution in terms of normal modes: 



Fi{p)e 



fi{0,p,z) 



i{e- 



Flip)' 



-i{e-LO*z) 



Al e + Xi e 



iY* 



(A5) 
(A6) 
(A7) 
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where the symbol * refers to the complex conjugate and 
in general to £ C. Making use of the above ansatz in 
the linearized system of equations returns the following 
consistency equation: 

uj^ r^dp^^ (A8) 

often referred to as to the dispersion relation. To de- 
termine whether a given distribution /o(p) is stable or 
unstable, one can solve the above dispersion relation and 
estimate the sign of the imaginary part of w. Depend- 
ing on the sign the field grows exponentially (instability) 
or oscillates indefinitely (stability). If the selected initial 
condition is parametrized via an adjustable parameter, 
one can then calculate the corresponding theshold value 
which discriminates between stable and unstable regimes. 
This analytical procedure can be persecuted in simple 
cases, as the one addressed in this paper (the watcrbag). 
For more complicated situations one can resort to the 
celebrated Nyquist method, first introduced in plasma 



physics [13 (see also [l^). For the case at hand, /o(p) 
takes the form: 



h{p) = 2^ 2A^ + ^P) - ^(P - ^P)] (A9) 

where 9 stands for the Heaviside function. Inserting (|A9[) 
into (|A8|) . carrying out the integral explicitly and looking 
for the value of Ap which sets the transition between 
complex and real w, leads to Ap ~ 1.37 or equivalently 
e = (Ap)V6 ~ 0.315. 
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